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Abstract 



In the framework of coupled-channel formalism the relativistic four-quark equations are 
found. The dynamical mixing of the meson-meson states with the four-quark states is 
considered. The four-quark amplitudes of the tetraquarks with open charm, including 
u, d, s, c quarks, are constructed. The poles of these amplitudes determine the masses 
and widths of tetraquarks. 
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I. Introduction. 

The discovery by the Belle Collaboration [1] of the very narrow X(3872) resonance 
has triggered the interest in the charmonium-like states, both experimentally and theo- 
retically. The observation of X(3872) has been confirmed by CDF [2], DO [3] and BaBar 
Collaboration [4]. Belle Collaboration observed the X(3940) in double- charmonium 
production in the reaction e + e" — > J/ip + X [5]. The state, designated as X(4160), 
was reported by the Belle Collaboration in Ref. 6. The fact that the newly found states 
do not fit quark model calculations [7] has opened the discussion about the structure 
of such states. Maiani et al. advocate a tetraquark explanation for the X(3872) [8, 
9]. On the other hand, the mass of X(3872) is very close to the threshold of D*D 
and, therefore, it can be interpreted as molecular state [10 - 15]. Several review pa- 
pers, as for example [16, 17], discuss the difficulty of interpreting these resonances as 
charmonium states. 

In the present paper the relativistic four-quark equations are found in the framework 
of coupled-channel formalism. The dynamical mixing between the meson-meson states 
and the four-quark states is considered [18 - 20]. Taking the X(3872) and X(3940) as 
input [21] we predict the masses and widths of S'-wave tetraquarks with open charm 
(Table I). 

II. Four-Quark Amplitudes for the Tetraquarks with Open Charm. 

We derive the relativistic four-quark equations in the framework of the dispersion 
relations technique. 

The correct equations for the amplitude are obtained by taking into account all 
possible subamplitudes. It corresponds to the division of complete system into subsys- 



tems with the smaller number of particles. Then one should represent a four-particle 
amplitude as a sum of six subamplitudes: 



A = A 12 + A 13 + A U + A 23 + A 24 + A U . (1) 

This defines the division of the diagrams into groups according to the certain pair 
interaction of particles. The total amplitude can be represented graphically as a sum 
of diagrams. 

We need to consider only one group of diagrams and the amplitude corresponding, 
for example A 12 . The relativistic generalization of the Faddeev-Yakubovsky approach 
[22, 23] for the tetraquark is obtained. We shall construct the four-quark amplitude of 
cuuu tetraquark in which the quark amplitudes with quantum numbers of 0~ + and 1" 
mesons are included. The set of diagrams associated with the amplitude A 12 can fur- 
ther be broken down into four groups corresponding to subamplitudes: Ai(s, s± 2 , s 34 ), 
A 2 (s, s 23 , su), A 3 (s, s 23 , S123), At(s, S14, S124), if we consider the tetraquark with the 
spin-parity J pc = ++ (cuuu). 

Here is the two-particle subenergy squared, corresponds to the energy 
squared of particles i, j, k and s is the system total energy squared. 

In order to represent the subamplitudes A^s, s 12 , s 34 ), A 2 (s, s 23 , s 14 ), A 3 (s, s 23 , s 123 ), 
A 4 (s, S14, S124) in the form of dispersion relations it is necessary to define the ampli- 
tudes of quark- ant iquark interaction a n (sik). The pair quarks amplitudes qq — > qq are 
calculated in the framework of the dispersion N/ D method with the input four-fermion 
interaction [24 - 26] and the quantum numbers of the gluon [27]. The regularization 
of the dispersion integral for the D-function is carried out with the cutoff parameter 
A. The four-quark interaction is considered as an input [27]: 

9v (qMfl^q) 2 + g { v (qMf-f^q) (sX-f^s) + g^ s) (sX-f^s) 2 . (2) 

Here If is the unity matrix in the flavor space (u,d). X are the color Gell-Mann 
matrices. Dimensional constants of the four-fermion interaction g v , and gy^ are 
parameters of the model. At gy = gft = gy^ the flavor SU (3) f symmetry occurs. The 
strange quark violates the flavor SU(3)f symmetry. In order to avoid an additional 
violation parameters, we introduce the scale shift of the dimensional parameters [27]: 

9 = ^=^2^=^. (3) 

a - (4) 

(mi + m k y 

Here mj and m k are the quark masses in the intermediate state of the quark loop. 
Dimensionless parameters g and A are supposed to be constants which are independent 
of the quark interaction type. The applicability of Eq. (2) is verified by the success of 
De Rujula-Georgi-Glashow quark model [28], where only the short-range part of Breit 
potential connected with the gluon exchange is responsible for the mass splitting in 
hadron multiplets. 

We use the results of our relativistic quark model [27] and write down the pair 
quarks amplitude in the form: 



(5) 



(m i +m k ) 2 h 

B n (s tk )= / j^ ^W^W . ( 6 ) 

(mi+m k y 

Here G n (si k ) are the quark- ant iquark vertex functions. The vertex functions are de- 
termined by the contribution of the crossing channels. The vertex functions satisfy 
the Fierz relations. All of these vertex functions are generated from gy, gy' 1 and 
gy S \ B n (s ik ), p n (s ik ) are the Chew-Mandelstam functions with cutoff A and the phase 
spaces, respectively. 

Here n — 1 determines a gg-pairs with J pc = 0~ + in the l c color state, n = 2 
corresponds to a gg-pairs with J pc = 1 in the l c color state, and n = 3 defines the 
gg-pairs corresponding to tetraquarks with quantum numbers: J pc = ++ , 1 ++ , 2 ++ . 

In the case in question, the interacting quarks do not produce a bound state; there- 
fore, the integration in Eqs. (7) - (10) is carried out from the threshold (m, + m k ) 2 
to the cutoff A(i, k). The coupled integral equation systems (the tetraquark state with 
n = 3 and J pc = ++ for the cuuu) can be described as: 

A i \ ^iB 2 (si 2 )B 2 (s u ) . i / O <~)\ A f i i \ 

Al(s,S 12 ,S U ) = 77 5-7 7777 777 77 +2J 2 (S 12 ,S34,2,2)A 3 (S,S 2 3,S 12 3) 

+ 2J 2 (S12,S34,2,2)A 4 (S, s'i4,s'i 24 ) , (7) 

A I \ ^2-Bi(s 23 )5i(Si 4 ) „J, ' ' \ 

^i2lS, S23, Sl4j - [1 _ _B 1 ( S23 )][1 _ 5 1 (s 14 )] ^ J2 ^ 23 ' Sl4 > X > L ) /i 3{S, S 34 , S 234 J 

+ 2J 2 (s 23 ,s 14 , l,l)A 4 (s,s 12 ,s 123 ) , (8) 



A 3 (s, s 23 , si 23 ) = A3jg 3(g2s) + 2 j 3 ( S235 3)Ai(s, s' 12 , s' u ) + J 3 (s 23 , 3)A 2 (s, s' 12 , s 



+ Jl( S 23: 3)A 4 (S, S 34 , S 2 34) + -A(«23, 3)A 3 (s, 

S 12' S 123 ), (9) 

A 4 (s, su, s i24 ) = J^~^T^\\ + 2 ^(si4, 3)Ai(s, s' 13 , s 24 ) + 2J 3 (si4, 3)A 2 (s, s' 13 , s 24 ) 

[1 - ±* 3 (Sl 4 JJ 

S 14> S 134 

) + 2J 1 (s 14 ,3)A 4 (s, 

S 14) S 134 ) , (io) 

where Aj, i — 1,2,3,4 are the current constants. They do not affect the mass spectrum 
of tetraquarks. We introduce the integral operators: 



f(. i\ - G i( s u) [ ds u Gi{s 12 )pi{s 12 ) tdz x ( . 

Ji\s 12 ,i) - 7- —, — 77 / ; / — , {i-i- 

[1-Bi(s 12 )\ J 7T S 12 -Si 2 J 2 
(mi+m 2 ) 2 



j ( j x Gi(s 12 )G p (sm) f ds' 12 Gi(s' 12 )pi{s' 12 ) 



[1 - Bi(s 12 )][l -B p (s u )} s' 12 -s 12 



(m3+m 4 ) 2 A 

4 , , _ , , s / , s +1 , +1 



X 



r ds' u G p (s' u )p p (s' u ) r dz 3 r dz 4 
J 7T s' 34 - s 34 / 2 i 2 ' 



^3(^12,0 



(m 3 +m,4) 2 

G,(si 2 ,A) I / ds? 12 G l (s' 12 ,A)p l (s' 12 ) 



(m 1 +m 2 ) 2 A 



X 



[1 - S,(S12, A)] 47T 7 7T s' 12 - S12 

(mi+m 2 ) 2 

+1^ +1 4 ^ 

/ dz dz 2 ; , (13) 

i 2 i /_ v/l - ^ 2 - ^ - 5 + 2^ 



here /, p are equal to 1 — 3. 

In Eqs. (11) and (13) Z\ is the cosine of the angle between the relative momentum 
of the particles 1 and 2 in the intermediate state and the momentum of the particle 3 
in the final state, taken in the cm. of particles 1 and 2. In Eq. (13) z is the cosine 
of the angle between the momenta of particles 3 and 4 in the final state, taken in the 
cm. of particles 1 and 2. z 2 is the cosine of the angle between the relative momentum 
of particles 1 and 2 in the intermediate state and the momentum of the particle 4 in 
the final state, is taken in the cm. of particles 1 and 2. In Eq. (10): z 3 is the cosine 
of the angle between relative momentum of particles 1 and 2 in the intermediate state 
and the relative momentum of particles 3 and 4 in the intermediate state, taken in the 
cm. of particles 1 and 2. z 4 is the cosine of the angle between the relative momentum 
of the particles 3 and 4 in the intermediate state and that of the momentum of the 
particle 1 in the intermediate state, taken in the cm. of particles 3, 4. 

We can pass from the integration over the cosines of the angles to the integration 
over the subenergies [29]. 

Let us extract two-particle singularities in the amplitudes Ai(s, s± 2 , S34), A 2 (s, s 23 , su), 
A 3 (s, S23, S123) and A 4 (s, s 14 , Si 24 ): 

At \ Oii(s, s ik , sim)B 2 (s ik )B 2 (si m ) 

A 1 (s,s i k,s lm ) = - r , w —, r r , (14) 

[1 - B 2 (s ik )\[l - B 2 {s lm )\ 

A( s _ a2(s,s i k,s lm )B 1 (s i k)B 1 (si m ) 

A 2 {s, s ik , si m ) — — — — — — — , [lb) 

[1 - Bi(s ifc )J[l - Si(s Jm )J 

A j (s,s ik ,s m ) = — — r , j = 3-4. (16) 

1 - B 3 {s ik ) 

We do not extract three-particles singularities, because they are weaker than two- 
particle singularities. 

We used the classification of singularities, which was proposed in paper [30]. The 
construction of the approximate solution of Eqs. (7) - (10) is based on the extraction of 
the leading singularities of the amplitudes. The main singularities in Si k ~ (m; + m k ) 2 
are from pair rescattering of the particles i and k. First of all there are threshold square- 
root singularities. Also possible are pole singularities which correspond to the bound 
states. The amplitudes apart from two-particle singularities have triangular singulari- 
ties and the singularities defining the interactions of four particles. Such classification 
allows us to search the corresponding solution of Eqs. (7) - (10) by taking into account 
some definite number of leading singularities and neglecting all the weaker ones. We 



consider the approximation which defines two-particle, triangle and four-particle singu- 
larities. The functions a^s, s 12 , s 34 ), a 2 (s, s 23 , s 14 ), a 3 (s, s 23 , s 123 ) and a 4 (s, s 14 , s 124 ) 
are the smooth functions of Sik, Siki, s as compared with the singular part of the ampli- 
tude, hence they can be expanded in a series in the singulary point and only the first 
term of this series should be employed further. Using this classification, one defines 
the reduced amplitudes a±, a 2 , a 3 , a 4 as well as the £?-functions in the middle point 
of physical region of Dalitz-plot at the point s : 



Sq = 0.25(771; + m k y.s , 



s + 2 £ m 2 



s 123 = 0.25s E (mi + m k ) 2 -^ra, 2 , s = j— ^ . (17) 

%f i=l 0.25 E {m + m k f 

i,k=l 

Such a choice of point s allows us to replace integral Eqs. (7) - (10) by the algebraic 
equations (18) - (21) respectively: 

a ± = Ai + 2a 3 JB 1 (2, 2, 3) + 2a 4 JB 2 (2, 2, 3) , (18) 
a 2 = \ 2 + 2a 3 JB 3 (l, 1, 3) + 2a 4 JB 4 (l, 1, 3) , (19) 
a 3 = X 3 + 2«i Jd(3, 2, 2) + 2a 2 JC 2 (3, 1, 1) + a 4 JAi(3) + a 3 JA 2 (3) , (20) 

a 4 = A 4 + 2«! JC 3 (3, 2, 2) + 2a 2 JC 4 (3, 1, 1) + a 3 JA 3 (3) + a 4 JA 4 (3) . (21) 

We use the functions JA^l), JBi(l,p,r), JCi(l,p,r) (l,p,r = 1 — 3), which are 
determined by the various s 1 q (Eq. 17). These functions are similar to the functions: 

G 2 {sf)B 2 {sf) ~^~ds' 12 Pl (s' 12 ) + rd Zl 1 



JMl) " A(4 2 ) , J „ tt s ; 2 - Sl J 2 i:22i 



(mi+m2) 2 — 1 



>23y 



, R r , , g^g|mg^ > ds' 12 Pl { s ' 12 ) 



B^B^sf) J 7T s' 12 -s 12 

(mi+m 2 )^ 



(mQ+m 4 ) 2 A 



r ds u Pp {s u ) r dz 3 r dz 4 1 
7 7T s' 34 -s 34 J 2 i 2 1-B r (4)' 

(m 3 +m 4 )- i! -1 -1 



, fi 2 (^ 2 ,A)73 p (^)73 r (^)l-^(4 2 ) 1 f ds^_Pi{s' 12 ) 
3[ ,P ' ] ~ l-B,(s\ 2 .A) B,(sn 2 ) An J n s\ 



1 - Bi(sq 2 ,A) B^sl 2 ) 4tt J n s' 12 -s 12 

(m\+m2) 



+1 , +1 2 2 

X 



z 2 -\- 2zz\z 2 



1 



(24) 



x 



[1 - B p (s> 23 m ~ B r (s' u )} ' 



A(ik) 



A(ik), 



if A(ik) < (v^i + m 3 ) 2 



(25) 



(y/SV23 +m 3 ) 2 , 



if A(ik) > (^3~ + m 3 ) 2 



The other choices of point s do not change essentially the contributions of a\, a 2 , 
a 3 and a±, therefore we omit the indices sff. Since the vertex functions depend only 
slightly on energy it is possible to treat them as constants in our approximation. 

The solutions of the system of equations are considered as: 



where zeros of D(s) determinants define the masses of bound states of tetraquarks. 
Fi(s, A,) determine the contributions of subamplitudes to the tetraquark amplitude. 

III. Calculation results. 

Our calculations do not include the new parameters. We use the cutoff A = 10.0 and 
the gluon coupling constant g = 0.794, which are determined by fixing the tetraquark 
masses for the states with the hidden charm J pc = 1++ X(3872) and J pc = 2 ++ X(3940) 
[21]. The quark masses of model m u>d = 385 MeV and m s = 510 MeV coincide with 
the ordinary meson model ones [27]. In order to fix anyhow m c = 1586 MeV, we use 
the tetraquark mass for the J pc = 2 ++ X(3940). The masses and widths of meson- 
meson states with the spin-parity J pc = ++ , 1 ++ , 2 ++ are given in Table I. In our 
paper we predicted the scalar tetraquark with the mass M = 2610 MeV and the width 
T ++ = 180 MeV (channel D°i]). We calculated the scalar tetraquark with the mass 
M = 2691 MeV and the width r ++ = 110 MeV (channels D+r] and D°K + ). The 
other scalar tetraquark is predicted as M = 2805 MeV and the width T ++ = 80 MeV 
(channel D+K~). The tetraquarks with the spin-parity J pc = 1 ++ , 2 ++ (Table I) have 
only the weak decays. 

The functions F^s, Aj) (Eq. (26)) allow us to obtain the overlap factors / for 
the tetraquarks. We calculated the overlap factors / and the phase spaces p for the 
reactions X — > M X M 2 (Table II). The widths of the scalar tetraquarks with open 
charm are obtained (Table I). We considered the formula r ~ f 2 x p [31], there p is 
the phase space. The widths of the tetraquarks are fitted by the fixing width T 2 ++ = 
(39 ± 26) MeV [21] for the tetraquark X(3940) (ccuu) with the spin-parity J pc = 2++. 

In the open charm sector the scalar tetraquarks have relatively small width ~ 
100 — 200 MeV, so in principle these exotic states could be observed. The low-lying 
tetraquarks with the open charm were calculated in other works [32, 33]. These states 
appear as narrow states. In our model the tetraquarks with the open charm and the 
spin-parity J pc = 1 ++ , 2 ++ can decay only in the weak channels. 
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Table I. Masses and widths of tetraquark with open charm. 

Parameters of model [21]: quark masses m u ^ = 385 MeV, m s = 510 Me V and m c = 
1586 MeV; cutoff parameter A = 10.0, gluon coupling constant g = 0.794. 



Tetraquark 


j V c 


Mass {MeV) 


Width {MeV) 


j V c 


Mass (MeV) 


jpc 


Mass {MeV) 


(cu)(uu) 


0++ 


2610 


180 


1++ 


2672 


2++ 


2736 


(cs)(uu) 
(cu)(us) 


0++ 


2691 


110 


1++ 


2770 


2++ 


2851 


(cu)(ss) 
(cs)(su) 


0++ 


2805 


80 


1++ 


2890 


2++ 


2975 



Table II. Overlap factors / and phase spaces p of tetraquarks with open charm. 



Tetraquark 


(channels) 


jpc 


/ 


P 


(cu)(uu) 




0++ 


0.396 


0.325 


{cs){uu) 
{cu){us) 


D+r, 
D°K+ 


0++ 


0.246 
0.183 


0.300 
0.414 


(cu)(ss) 
{cs){su) 


D% 
DtK- 


0++ 


0.192 
0.237 


0.407 
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